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^ ' Abstract. We demonstrate a technique for exciting spin waves in an ultracold gas of 

^r\. ^^Rb atoms based on tunable AC Stark potentials. This technique allows us to excite 

normal modes of spin waves with arbitrary amplitudes in the trapped gas, including 
C ■ dipole, quadrupole, octupole, and hexadecapole modes. These modes exhibit strong 

2 , nonlinearities, which manifest as amplitude dependence of the excitation frequencies 

qH' and departure from sinusoidal behavior. Our results are in good agreement with a full 

treatment of a quantum Boltzmann transport equation. 
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Spin waves in ultracold atomic gases offer a new regime for exploring quantum dynamics. 



G> ! Strikingly, applying a brief coupling between two internal quantum states in an otherwise 

^r^ , classical gas can excite macroscopic collective quantum dynamics that persist for orders 

"^ ! of magnitude longer than the initial perturbation. Furthermore, many experiments 

L^ ! must go to great lengths to achieve nonlinear couplings for studying more complicated 

O ! behaviors than those present in simpler linear systems, whereas in the ultracold gas 

system, it is nearly impossible not to drive the system into a nonlinear regime. 

The first observations of spin waves in a fluid were in spin-polarized H and later 
rS I in liquid ^He and ^He-^He mixtures [1]. The high density and strong coupling between 

5^ i particles meant that these systems were deep in the hydrodynamic regime of the spin 

wave spectrum. Thus, they gave rise only to linear excitations, and the spin orientations 
never diverged far from equilibrium. The discovery of spin waves in ultracold gases 
opened the door to more detailed studies of these excitations due to the low density, 
ease of manipulation, and imaging possibilities of such gases [21 13]. It was immediately 
obvious that spin waves in ultracold gases were almost intrinsically nonlinear, since 
they occurred between the hydrodynamic and coUisionless limits. They showed large 
deviations from spin equilibrium, the spin states almost fully separating spatially, 
and strong amplitude dependence on the spin wave frequencies. These nonlinearities 
complicated analyses, as simplifying linear approximations were not appropriate. Much 
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care was taken to eliminate the effect of nonlinearities by extrapolating to low amplitude 
limits, but the nonlinear behavior itself was not studied. 

In this work, we specifically study the nonlinear nature of spin waves in an ultracold 
trapped alkali gas. We do so by controllably exciting different spin wave modes with 
arbitrary amphtude and studying the frequency shifts produced. Previous work on spin 
waves in trapped atomic gases relied on parameters intrinsic to the trapped atom system 
to excite spin waves - in particular density- dependent mean-field shifts as well as the 
optical or magnetic trapping potential itself [SI HJ |5] . Here we eliminate this driving 
mechanism and instead add an external optical potential relying on the AC Stark shift, 
so that spin wave amplitudes are decoupled from the trap parameters and arbitrary 
modes can be excited. 

With this new method, we not only excite the quadrupole mode (mode number 
d = 2) observed in [3], but we also demonstrate the flexibility of optical excitation 
of spin waves by driving dipole, octupole, and hexadecapole modes {d = 1,3, and 4 
respectively), which have not been previously observed in ultracold gas systems. We 
excite these modes with varying amplitudes and show that the excitation frequencies 
depend strongly on the amplitude, so much so that in some instances it is virtually 
impossible to excite spin waves in the linear regime where there is no amplitude 
dependence. This behavior is in good agreement with numerical simulations of a 
Boltzmann transport equation. Finally, we highlight the nonlinearity of this system 
by demonstrating that even after short time periods after beginning the excitation, 
during which the spin perturbations have not yet become large, the spin dynamics are 
inconsistent with those predicted by a simplified linear analysis. 

2. Spin wave theory 

Spin waves are macroscopic, collective excitations of atomic spin vectors that arise 
from the identical spin rotation effect (ISRE) due to exchange scattering between 
indistinguishable particles in an inhomogeneous potential ^. The spin system we 
consider is a harmonically confined ultracold gas of spin-1/2 atoms, with spin states 
|1) and |2) coupled electromagnetically. Atoms in an equal coherent superposition of 
the two spin states, \-ip) = (|1) + e*'^|2))/A/2, in a spatially varying differential potential 
Udis{f) will acquire a spatially dependent phase, 0(r, t) = Udm{r)t/h. This phase 
represents a transverse rotation of the local spin vector S{r). As an atom oscillates 
in the confining potential, it undergoes coherent exchange scattering with other atoms 
having slightly different spin vectors. Indistinguishability of quantum particles requires 
the scattered wavefunction to be symmetrized with respect to forward (unscattered) 
and backwards scattered terms, which leads to a rotation of the spin vectors of the two 
interacting particles about their combined total spin. This ISRE has been observed in 
both bosonic and fermionic systems OH]. The ISRE gives rise to coherent spin currents, 
and the atomic spins undergo spatiotemporal oscillations of the spin vectors, i. e. spin 
waves, about their equilibrium spin configuration as determined by f/diff- 
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Spin dynamics in this system are described by a Boltzmann spin transport equation 
[TJ m |9]. Since the geometry of our atom trap is highly elongated in one direction (z), 
we consider only the one- dimensional evolution of the spin position and momentum 
distribution a{p,z,t) and the experimentally accessible total spatial spin distribution 
S{z,t) = J dp a{p,z,t)/2nh. The transverse components of S, Su and S'^, represent 
the spin coherence or the superposition phase 0, and the longitudinal spin component 
represents the spin differential between states |1) and |2), 5*^; = A^i — N2. We adopt the 
notation of ^U\ and write the quantum Boltzmann equation as 



dcr p dcr dUcxt dcr dcr 

dt m dz dz dp dt 



(1) 



ID 

where we have suppressed the explicit dependence on (p, z, t) for clarity. The first 
three terms are the total time derivative of cr, Ue-xt is the harmonic magnetic trapping 
potential, and the right hand term represents a radially averaged 1-D collision term 
proportional to the elastic scattering rate. We employ the analytical relaxation time 
approximation for the collision integral derived in [10]. The remaining term is a spin 
torque that couples the components of the spin vector via the differential potential and 
mean field coupling, and is given by 

n = {U^^w^ + gS)/n, (2) 

with g = Anh^a/m for mass m and s-wave scattering length a|j It is this term that 
drives the spin rotations and provides the nonlinearity. 

To make the physics behind these spin oscillations more accessible, ^U\ performed 
a small amplitude moment expansion of the spin distribution, |7] considered the 
excitations in the hydrodynamic limit where spin perturbations are never far from 
equilibrium values, and pJJ compared the two methods. These linear approximations 
give analytical results that allow for a better conceptual understanding of the 
phenomenon, and they accurately predict the small amplitude behavior of the spin 
excitations. However, linear approximations fail to capture the full dynamics of large 
amplitude spin waves, and equation ([T]) must be treated numerically as was done in 
[H [9]. The observed nonlinear spin wave behaviors that are predicted by numerical 
solutions to equation ([T]) include amplitude dependent frequencies, frequency chirps as 
the excitations damp to equilibrium, and nonsinusoidal initial growth of S^ at the onset 
of the excitations. 

3. Spin wave apparatus 

The system in which we study nonlinear spin waves consists of two hyperfine ground 
states of ^'^Rb, which comprise a pseudo-spin 1/2 doublet. The two states, |1) = |-F = 

I The scattering lengths a^ between states \i) and |j) (i,j = 1 or 2) vary by only a few percent in 
^^Rb. This difference has only a minor effect on the dynamics and is not important for effects described 
herein, in contrast to phase separation in Bose-Einstein condensates, which occurs to minimize mean- 
field energy. 
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Figure 1. Schematic of tlie apparatus used for optical excitation of spin waves (not to 
scale), including an absorption image of the |1) state near the peak of a hexadecapole 
oscillation featuring three local maxima [S rotated towards |1)) and two local minima 
{S rotated towards |2)). 



l,mp = —1) and |2) = \F = 2,mp = 1), are coupled via a two-photon microwave 
transition using a 6.8 GHz microwave photon and a 3.14 MHz radiofrequency (rf) 
photon, both detuned 700 kHz from the intermediate \F = 2, Trip = 0) state. The 
two-photon Rabi frequency is 3.2 kHz. The atoms are confined in a hybrid loffe- 
Pritchard magnetic trap similar to that described in [12]. The trap is cylindrically 
symmetric and elongated in the axial (z) direction with an aspect ratio of 37:1 (axial and 
radial frequencies 6.7 Hz and 247 Hz respectively), which allows us to radially average 
the spin profiles and treat the dynamics as one-dimensional. We use rf evaporative 
cooling to prepare a sample of atoms at temperature T = 650 nK and peak density 
uq = 2.5 X 10^^ m~^. These conditions are ~ 1.5 times the critical temperature for Bose- 
Einstein condensation, and the external degrees of freedom are well described as an ideal 
gas if one considers the gas as a single species and ignores the spin degree of freedom. 
The spin distribution for all experiments described herein is initialized by placing the 
atoms into a coherent superposition of |1) and |2) via a 7r/2 pulse, corresponding to a 
spatially uniform spin state S{z) = n{z)u. 

The differential potential experienced by the atoms has three contributions, as 
noted above: the magnetic confining potential, the mean field shift, and an additional 
optical potential. By tuning the magnetic field, the differential mean field shift can be 
almost exactly canceled by the Zeeman shift to give a fiat relative potential (see [2]). In 
the absence of any other inhomogeneities, no spin waves will arise in this configuration 
(although spin waves can still propagate in a uniform potential). 

We add a focused, spatially modulated laser to create tunable differential optical 
potentials and excite arbitrary spin wave modes (figure [1]). The AC Stark shift for a 
far detuned laser creates an optical dipole potential [/dip ~ I{r)/A, for laser intensity 
I{r) and detuning A. We use a 120 mW diode laser tuned A ^ 0.3 nm above the 
^^Rb cooling transition at 780.2 nm, focused to a beam waist of 40 /im. The 6.8 GHz 
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detuning difference between spin states is sufficient to produce a differential potential 
up to Udis/h ~ 500 Hz. A was chosen to minimize spontaneous emission but still keep 
it close enough to resonance to obtain the desired differential potentials without using 
high intensities, which could exert significant optical forces on the atoms. We verify this 
condition by noting that when all of the atoms are in a single spin state, the optical 
potential produces no noticeable change in the density distribution over the time scale of 
our experiment. Because of the elongated aspect ratio of the trap, all of the observable 
spin dynamics occur only in the axial direction, and the dipole potential can be radially 
averaged to produce a one-dimensional distribution. Thus, we need only one axis of 
spatial control, which we obtain by deflecting the laser with an acousto-optic modulator 
(AOM) as shown in figure [TJ 

We tailor spatially varying potentials with arbitrary symmetry by frequency and 
amplitude modulating (FM and AM) a 120 MHz rf signal applied to the AOM. 
The modulation rate of several kHz is significantly faster than all other time scales, 
including the trap frequencies and collision rates, and this potential can be viewed as 
a time-averaged, static optical potential. We create potentials with dipole, quadrupole, 
octupole, and hexadecapole symmetries {d = 1, 2, 3, and 4 respectively) by adjusting 
the AM and FM frequencies and phases relative to each other. For example, to create a 
dipole differential potential, we use 4 kHz phase- synchronized AM and FM, so that when 
the dipole beam is swept to one side of the cloud, its amplitude is maximal, while its 
amplitude is a minimum at the opposite end of the frequency/spatial sweep. Changing 
the modulation frequency of the AM and FM relative to each other while also shifting 
their relative phase produces higher d potentials. The overall amplitude of the potential 
is controlled via the laser power. Figure |2]^a)-(d) shows examples of potentials created 
in this manner. 

4. Optical excitation of spin waves 

To study the intrinsic spin wave modes of the system, we optically imprint these 
transverse spin profiles (phase profiles) onto the spatially uniform coherent superposition 
immediately following the 7r/2 pulse. We apply a 5 ms phase imprinting pulse with 
the same symmetry as the multipole mode we wish to drive, before turning off the 
optical potential and allowing the spin wave to oscillate freely in an unperturbed uniform 
potential. After a variable delay time, we measure the 1-D spatial distribution of the 
longitudinal spin vector S{z) by imaging the spatial distribution of the |1) and |2) 
states. We divide the cloud axially into equally sized bins and radially average the 
atomic distribution in each bin. Figures |2]^e)-(h) show spatiotemporal oscillations of Sy^ 
for the first four spin wave modes, plotted as a function of axial position in the cloud 
and free evolution time. The longitudinal spin is normalized by the population of each 
bin, Sw{z)/n{z) = {Ni — N2) / {N2 + Ni), to remove number fluctuations and account 
for trap loss. 

In all cases, the initially transverse spin vectors rapidly develop fluctuations of 
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Figure 2. L/(jiff/27rft for the (a) dipole, (b) quadrupole, (c) octupole, and (d) 
hexadecapole modes. C/diff was measured using a 7r/2 — 7r/2 Ramsey interferometer 
by varying the 7r/2 pulse separation for times up to 5 ms (before atoms will have 
moved appreciably in the axial direction), measuring the number of atoms returning 
to the |1) state, dividing the cloud into axial bins for radial averaging, and fitting each 
bin to a sinusoid to extract the local differential potential energy. Solid lines are to 
guide the eye. Each d-pole potential has d nodes and is centered and symmetric or 
antisymmetric, (e)-(h) show false color images of the spatiotemporal evolution of the 
normalized longitudinal spin distribution Sw/n for the above multipole modes. Data is 
interpolated between bins and time steps. Red/blue regions represent spatial locations 
where S^ has rotated out of the transverse plane and represent spin-state segregation. 



the longitudinal spin component, which manifests as spin segregation [2]. Higher order 
modes show proportionally higher frequencies due to the symmetry of the differential 
potentials. The maximum frequency for any mode d is governed by the trapping 
frequency according to /d = duz/^ii, which occurs in the zero density, collisionless 
limit. At higher density, /d is still proportional to d but also scales inversely with the 
exchange collision frequency, which is proportional to density. 

The damping rate, Fd, also increases with mode number d, as collisions bring to 
equilibrium faster those spin distributions that have greater inhomogeneity. These two 
differing timescales cancel out somewhat, and the quality factor, Qd = Syr/d/Fd, varies 
only slightly for the different modes. For the d = 2,3, and 4 modes, Qd ~ 1 ~ 4 and 
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Figure 3. Typical normalized spin wave moments, {Sw{z)z'^)n, for strongly driven 
excitations of the (a) dipole (• ) and quadrupole (■) modes, and (b) octupole (■) 
and hexadecapole (• ) modes. Solid lines are least-squares fits to damped sinusoids, 
for extracting frequencies and damping rates. 



the excitation damps after only one full oscillation. Because the dipole mode frequency 
is so low, /i ~ 1 Hz, it damps via an additional mechanism - dipolar relaxation of 
the 1 2) atoms leading to trap loss and additional shrinking of the spin vector - and 
Qi is smaller. It is important to note that while we distinguish between the different 
normal modes of the system based on their symmetry, these modes are not closed and 
can couple to other spin wave modes. This coupling is a significant contribution to the 
damping rate and is a further indicator of the nonlinearity of the system. 



4.1. Amplitude dependence of frequencies 

To highlight the nonlinear nature of these large amplitude spin waves, we systematically 
vary the amplitude of the initial perturbation to the transverse spin by adjusting the 
amplitude of the 5 ms phase imprinting pulse and studying the resulting spin wave 
frequencies (figure [3]). The multipole moment of each excitation is calculated according 
to {Sw{z,t)z'^)n to extract the oscillation frequencies, where the notation (. . .)„ denotes 
a density weighted average. This density weighting is necessary because the signal-to- 
noise ratio is lowest on the edges of the cloud where z'^ becomes large. Damped sine 
waves are then fit to the multipole moments to determine fd and r^;. Figure H] shows 
the dipole and quadrupole mode frequencies as a function of the amplitude of the initial 
perturbations. Negative values in figure H] are achieved by changing the phase of the FM 
and AM applied to the phase imprinting AOM. Other than reversing the sense of the 
rotation, equation ([T]) is unchanged by reversing the sign of dUdis/dz and d'^Udm/dz^. 

For all modes, fa increases significantly with amplitude, nearly doubling in the 
case of the quadrupole mode. As the frequency increases, however, F^ decreases faster. 
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Figure 4. Frequency of (a) dipole and (b) quadrupole spin wave oscillations as a 
function of excitation amplitude, parametrized by the gradient and curvature of U^^q 
at the center of the cloud. The solid lines are numerical simulations of the quantum 
Boltzmann equation (see text). Error bars are statistical and do not include systematic 
effects such as density calibration and the effects of nonsinusoidal oscillations. 



and Qd drops for larger perturbations. The octupole and hexadecapole modes are not 
studied in detail here, as they are quite sensitive to any imperfections in the optical 
driving potential. If the optical potentials do not have the exact symmetry of the 
desired mode, the d = 3 or 4 modes quickly damp and degenerate into longer lived 
(i = 1 or 2 excitations. 

To compare these measurements with the spin transport theory, we numerically 
simulate equation ([T]) using an alternating direction implicit finite difference method. 
We use the experimentally measured t/difr (see figure [2]^a) and (b) for instance) and 
account for imaging inefficiencies by adjusting the measured density by ~ 20% to match 
the low amplitude limit of the experimentally measured frequencies. For a given [/dig we 
calculate the multipole moment from the simulated S(^z^ t) and fit a damped sinusoid 
to it for comparison with the measured frequencies. As seen in figure HI the simulations 
show good agreement with the data, indicating that the evolution of the spin vector is 
well described by the Boltzmann equation. 

The primary systematic error in these measurements is loss of atoms from the 
1 2) state due to dipolar relaxation. As the spin states separate, dipolar relaxation is 
accelerated, significantly altering the dynamics. Furthermore, the driving potentials 
we use are not purely linear or purely quadratic, and lead to excitation amplitudes 
and patterns that are not completely characterized by a single parameter [dUdis/dz or 
d'^Udis/dz'^) taken at the center of the trap. It is important to note that the oscillations 
of Sw are not truly sinusoidal; however, since Qd ~ 3, a sinusoidal fit over a single 
oscillation is a reasonable characterization of the primary time scale. We also verify via 
Monte Carlo simulations that there is no systematic error from fitting damped sinusoids 
when /d and T^ are not well separated, so long as Q^ > 1. However, the fitted frequencies 
are systematically lowered when gamma becomes larger than than the frequency - this 
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effect is seen in figure Hl^b) in both the theory and data for the most strongly driven 
spin waves above 3 kHz/mm^ initial curvature. In this region the fits, while repeatable, 
are not true measures of an oscillation frequency. 

Despite the symmetry of the Boltzmann equation, quadrupole excitations in our 
experiment are asymmetric with respect to the sign of the curvature of t/difr, as seen 
in figure Hl^b). There are several reasons for this asymmetry. Although we parametrize 
the potentials by the curvature at the center of the trap, the quadrupole differential 
potential is not truly quadratic and rolls off near the edges of the cloud (figure El^b)); 
this roll off changes shape somewhat when we reverse the sign of t/dis- This effect is 
included in our simulations by using the full experimentally determined potentials as 
inputs, which accounts for some of the positive/negative asymmetry in the frequency. 
Another significant effect comes from, again, losses of the |2) state via dipolar relaxation 
collisions. When |2) segregates to the trap center (positive curvature), dipolar loss 
slightly speeds up the return of the |1) state to the center of the trap. However, for 
negative curvature, |2) segregates to the outside of the trap and preferentially leaves the 
trap, making the return oscillation appear slower since the spin current is not conserved. 
We account for this effect at each time step of the simulation by adding a small 
phenomenological loss to the projection of S onto |2) proportional to the local |2) state 
density. With these modifications, we obtain good agreement between the simulations 
and measurements, both showing strong nonlinear behavior. The linear regime, where 
the frequency curve flattens, occurs at extremely small excitation amplitude and is 
almost nonexistent on this scale. We are not able to drive small enough amplitude spin 
waves to explore this limit. 

4-2. Short time behavior 

The behavior of the spin vector shortly after the optical excitation commences also 
highlights the nonlinear nature of these spin waves. Both the small amplitude linearized 
moment method of [10] and the near equilibrium hydrodynamic approach of [7] predict 
that Sw initially grows as t^. However, equation ([1]) can be expanded for t < Tcoi, the 
mean collision time pi [13]. In this limit, the initial growth of Syj goes as t^: 

Sw gn ksT [d^Udis 2z dUdm 



n 48/i m 



t\ (3) 



dz"^ al dz 
where a^ is the axial Gaussian width of the atomic distribution. 

Figure |5] shows Sw in the center of the cloud for a hexadecapole excitation similar 
to that shown in figures WA) and (h) for times shorter than t^ox — 12 ms. For this 
study we use a strong optical excitation to maximize nonlinearities and measure 5*^^, at 
a short time t after the start of the optical excitation pulse, all of which occurs before 
free evolution time shown in figure [2l We bin pixels in the center 23 /xm of the cloud 
image for averaging and normalize S^ by the number of atoms contained within that 
bin. We fit both t^ and f* power laws to S^- Although the ~ 5% shot-to-shot number 
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Figure 5. (a) Nonlinear beliavior in tlie initial growth of Sw (normalized) at the center 
of the cloud. Nonlinear least-squares fits to quartic (solid line) and quadratic (dashed 
line) are shown, along with the prediction of equation ([3]), the Boltzmann transport 
equation (dotted line). All points are an average of three experimental cycles, and 
the first point shows typical statistical error bars, (b) Normalized residuals from the 
quadratic (• ) and quartic fits (■) with solid lines to guide the eye. The residuals show 
a clear systematic trend for the quadratic fit (reduced x^ — 3.8), while the quartic fit 
residuals (x^ — 1.3) indicate that a quartic increase in S^ occurs. 



fluctuations contribute significant noise to a small Syj, the initial spin growth clearly 
does not grow quadratically and is better described by a quartic rise. 

The interpretation for this behavior is that the superposition must first build up 
phase gradients before coherent spin currents can develop, and only then does the ISRE 
begins to convert transverse phase gradients into longitudinal spin rotations. This multi- 
step process takes longer to initiate but rapidly accelerates once sufficient phase gradient 
has accumulated in what [I3] likened to a "density shock" rather than hydrodynamic or 
collisionless flow. The observed initial growth rate is even faster than that predicted by 
equation ([3]), which is most likely explained by the fact that the spin vectors continue to 
evolve during the several hundreds of microseconds imaging sequence prior to the image 
acquisition itself, as well as possible small errors in the calibration of density and U^m 
used in the simulations. As spin vectors rotate significantly out of the transverse plane 
and Syj gets large, the rate of increase slows down and becomes quadratic, approaching 
the sinusoidal behavior seen in figure [3l 

5. Conclusion 



We have demonstrated a versatile technique for exciting arbitrary modes of spin waves 
with arbitrary amplitudes. We have used this optical excitation technique to highlight 
the nonlinear behavior of an ultracold trapped atomic gas. Because the typical operating 
regime for spin waves in this system tends to be squarely between the hydrodynamic and 
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collisionless limits, the behavior of spin waves is more comphcated than that described 
by hnear or small amplitude approximations. In this case, it is almost impossible to 
excite purely linear spin waves, as even modest perturbations drive the sensitive atomic 
system into large amplitude fluctuations with significant nonlinear effects, including 
amplitude and time dependent changes in excitation frequency, increased damping, 
coupling between modes, and deviations from pure sinusoidal behavior. Furthermore 
the nature of the ISRE is complicated enough that even in the collisionless limit, 
the momentum cannot always be integrated out of the Boltztmann equation, and the 
dynamics are still not simple |ll[5l[Tl]. In either case, the ISRE and spin waves continue 
to offer rich and interesting examples of out-of-equilibrium quantum systems. 
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